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Abstract 

The ATLAS experiment at the LHC has measured the centrality dependence 
of charged particle pseudorapidity distributions over \r)\ < 2 in lead-lead colli- 
sions at a nucleon-nucleon centre-of-mass energy of t/s^ = 2.76 TeV. In order 
to include particles with transverse momentum as low as 30 MeV, the data 
were recorded with the central solenoid magnet off. Charged particles were 
reconstructed with two algorithms (2-point "tracklets" and full tracks) using 
information from the pixel detector only. The lead-lead collision centrality was 
characterized by the total transverse energy in the forward calorimeter in the 
range 3.2 < \rj\ < 4.9. Measurements are presented of the per-event charged 
particle pseudorapidity distribution, dN^/drj, and the average charged parti- 
cle multiplicity in the pseudorapidity interval |?7| < 0.5 in several intervals of 
collision centrality. The results are compared to previous mid-rapidity mea- 
surements at the LHC and RHIC. The variation of the mid-rapidity charged 
particle yield per colliding nucleon pair with the number of participants is con- 
sistent with lower ^/s NN results. The shape of the dN^/drj distribution is found 
to be independent of centrality within the systematic uncertainties of the mea- 
surement. 



1. Introduction 

Collisions of lead (Pb) ions at the Large Hadron Collider provide an oppor- 
tunity to study strongly interacting matter at the highest temperatures ever 
created in the laboratory pQ. Measurements of the centrality dependence of 
charged particle multiplicities and of charged particle pseudorapidity densities 
in such ultra-relativistic nucleus-nucleus (A+A) collisions provide essential in- 
formation on the initial particle or entropy production and subsequent evolution 
in the created hot, dense matter. Results from the Relativistic Heavy Ion Col- 
lider (RHIC) over the centre-of-mass energy range from 19.6 to 200 GeV indicate 
that the multiplicity of charged particles per colliding nucleon pair has a mild 
dependence on the collision centrality and that the pseudorapidity dependence 
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of the charged particle yield near mid-rapidity is essentially centrality indepen- 
dent 0. The weak variation of the multiplicity per colliding nucleon pair with 
centrality at RHIC was initially found to be inconsistent with models such as 
HIJING [3j which includes a mixture of soft and hard scattering processes with a 
Pt cutoff on the hard scattering contribution at 2 GeV, or with a beam-energy- 
dependent cutoff in a more recent version 0]. In contrast, calculations based 
on parton saturation invoking fcx factorization were able to reproduce both the 
shape and centrality dependence of the RHIC charged particle pseudorapidity 
distributions [5j[6]. However, more recent theoretical studies indicate that Ut 
factorization may not be applicable to nucleus-nucleus collisions, and improved 
soft+hard models may be able to describe RHIC multiplicity measurements. 
At the same time, older hydrodynamical models (e.g. Ref [7j) have had some 
success describing the energy dependence of the total multiplicity as well as ra- 
pidity distributions of identified hadrons, although their domain of applicability 
is still not fully established. 

Detailed measurements of the centrality dependence of charged particle mul- 
tiplicities and pseudorapidity distributions at the LHC together with the earlier 
RHIC measurements could provide essential insight on the physics responsi- 
ble for bulk particle production in ultra-relativistic nuclear collisions. Because 
hard scattering rates increase rapidly with centrality and ^/s NN , the combined 
RHIC and LHC measurements should provide a strong constraint on the con- 
tribution of hard scattering processes to inclusive hadron production subject 
to uncertainties regarding the shadowing of nuclear parton distributions at low 
x. Measurements at the LHC can also provide a valuable test of recent parton 
saturation calculations that still claim to be able to describe inclusive parti- 
cle production in ultra-relativistic nuclear collisions [51 [3] . Previous measure- 
ments at the LHC |H] have already started addressing some of the physics 
raised above. In particular, those earlier measurements found a rapid rise in 
the particle multiplicity at the LHC compared to naive extrapolations of RHIC 
measurements and a variation of mid-rapidity charged particle multiplicity with 
centrality similar to that observed at RHIC. 

This paper presents the results of ATLAS [TU] measurements of the per-event 
charged particle pseudorapidity distribution, dN^/d-q, in ^/s NN = 2.76 TeV 
Pb+Pb collisions over \r)\ < 2 and as a function of collision centrality with the 
goal of testing and extending the results of the previous LHC measurements. 
In this paper, iV ch denotes the per-event number of charged primary particles 
measured in an interval of r), which is the particle pseudorapidity Q The mea- 
surement was performed with the solenoid off, thereby allowing detection of 
charged particles down to very low transverse momenta (pr ~ 30 MeV). 



1 ATLAS uses a right-handed coordinate system with its origin at the nominal interaction 
point (IP) in the centre of the detector and the z-axis along the beam pipe. The x-axis 
points from the IP to the centre of the LHC ring, and the y axis points upward. Cylindrical 
coordinates (r, <j>) are used in the transverse plane, </> being the azimuthal angle around the 
beam pipe. The pseudorapidity is defined in terms of the polar angle 9 as r) = — lntan(#/2). 
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2. Experimental setup and event selection 

The measurements presented here were obtained using the ATLAS inner 
detector which contains both silicon pixel and silicon strip detectors and 
the ATLAS forward calorimeters. The charged particle multiplicity is measured 
using the pixel detector [12] which consists of three layers of pixel staves in the 
barrel region, inclined at an angle of 20°, at radii of 50.5, 88.5, and 122.5 mm 
from the nominal beam axis. The typical pixel size is 50 /im x 400 /jm in 
4> — z, and an average occupancy of about 0.5% is observed for the innermost 
pixel layer in central Pb+Pb collisions. To limit low-px multiple scattering 
losses in detector material, the measurement has been restricted to the barrel 
portion of the pixel detector, corresponding to pseudorapidity values in the 
range \rj\ < 2. Collision vertex positions were obtained by full reconstruction 
of nominally straight charged particle trajectories in the pixel and silicon strip 
detectors followed by reconstruction of a single collision vertex from the full set 
of particle trajectories. To maintain uniform acceptance of the pixel detector 
for the multiplicity measurement the vertex was required to lie within 50 mm 
of the nominal centre of the ATLAS detector in the longitudinal direction. 

The data for the measurements presented here were collected with a 
minimum-bias trigger. This required a coincidence in either the two minimum- 
bias trigger scintillator (MBTS) detectors, located at ±3.56 m from the inter- 
action centre and covering 2.1 < |^| < 3.9, or two zero-degree calorimeters 
(ZDCs), located at ±140 m from the interaction centre and covering \r]\ > 8.3. 
The threshold on the analog energy sum in each ZDC was set below the sin- 
gle neutron peak. The offline analysis required the time difference between the 
two MBTS detectors to be |Ai| < 3 ns to eliminate upstream beam-gas inter- 
actions, a ZDC coincidence to efficiently reject photo- nuclear events [13], and 
a reconstructed vertex satisfying the selection described above. The measure- 
ments presented in this paper were obtained from a 10 hour data-taking run 
corresponding to an integrated luminosity of approximately 480 mb -1 . A total 
of 1631525 events passed the trigger, vertex, and offline selections. 

3. Centrality 

In heavy ion collisions, "centrality" reflects the overlap volume of the two 
colliding nuclei, controlled by the classical impact parameter. That overlap 
volume is closely related to the number of "participants" , the nucleons which 
scatter inelastically in each nuclear collision. While the number of participants, 
Ap ar t, cannot be measured for a single collision, previous studies at RHIC and 
the SPS have demonstrated that the multiplicity and transverse energy of the 
produced particles are strongly correlated with -/Vp art . Because of this, the aver- 
age number of participants can be accurately estimated from a selected fraction 
of the multiplicity or transverse energy distribution |14| . In ATLAS, the Pb±Pb 
collision centrality is measured using the summed transverse energy (EiJx) in 
the forward calorimeter (FCal) over the pseudorapidity range 3.2 < I77I < 4.9, 
calibrated at the electromagnetic energy scale. An analysis of the FCal E£/r 
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distribution after application of all trigger and selection requirements gives an 
estimate of the fraction of the sampled non-Coulomb inelastic cross section of 
/ = 98 ± 2%. This estimate was derived from comparisons of the measured 
FCal Y,Et distribution with a simulated Y,E^ distribution. The simulated dis- 
tribution was obtained from a convolution of ^/s = 2.76 TeV proton-proton data 
with a Monte Carlo (MC) Glauber calculation [HI [15] of the number of effective 
nucleon-nucleon collisions. This quantity was calculated as a linear combina- 
tion of the number of participants and the number of binary collisions, similar 
to what was done in a previous analysis [16] . The value of / and its uncertainty 
was estimated by systematically varying the effect of trigger and event selection 
inefficiencies as well as backgrounds in the most peripheral Ei?x interval. This 
was done by artificially injecting and removing counts in that interval in order to 
achieve the best agreement between the measured and simulated distributions. 
The estimate of / was made after removal of a 1% background contamination 
in the most peripheral events that was evaluated using comparisons of solenoid 
magnet-on and solenoid magnet-off data and which was attributed to photo- 
nuclear events. 

For the results presented in this paper, the minimum-bias FCal Si?x dis- 
tribution was divided into centrality intervals according to the following per- 
centiles: 10% intervals over 0-80%, 5% intervals over 20 — 80% and 2% intervals 
over — 20%. By convention, the 0-10% centrality interval refers to the 10% 
most central events - the events with the highest £Z?t values - and increasing 
percentiles refer to events with successively lower YiE^. The average number 
of participants, (iVp art ), was evaluated for each of the experimental centrality 
intervals by dividing the Glauber model ££/r distribution into the same per- 
centile centrality intervals used for the data and evaluating the average number 
of participants of the Glauber MC events contributing to a given interval. This 
procedure incorporates more realistic fluctuations into the estimation of (-/Vp art ) 
than would be achieved by binning in either -/V part itself or in the classical 
impact parameter. The systematic errors on (iV part ) were evaluated from the 
quoted uncertainty on / and the known uncertainties in the nuclear density 
parameters as well as the assumed total inelastic nucleon-nucleon cross section 
of ctnn = 64 ±5 mb Ql]. 

4. Reconstruction of charged particle multiplicity 

In the offline analysis, adjacent hits in the pixel modules were grouped into 
clusters using standard techniques. Two methods were, then, used to recon- 
struct charged particles from the pixel clusters. In one method, a Kalman 
Filter-based tracking algorithm, similar to that deployed in proton-proton col- 
lisions 18], was applied only to the pixel layers ("pixel tracks"). The other 
method, the "two-point tracklet" algorithm, used the reconstructed primary 
vertex and clusters on the first pixel layer to define a search region for clusters 
in the second layer consistent with a nominally straight track. Candidate track- 
lets were required to have deviations between projected and measured cluster 
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positions in the second pixel layer in pseudorapidity and azimuth, A77 and A<f>, 
respectively, satisfying 
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The widths of the Arj and Ac/) distributions characterized by the pseudorapidity- 
dependent resolutions <j v (tj) and o^if) were obtained from the MC simulations 
described below. The 77 and <j) values of the reconstructed tracklets were de- 
termined using the cluster position on the first layer and the primary vertex 
position. The two-point tracklet analysis excluded clusters with low energy de- 
posits inconsistent with minimum-ionizing particles originating at the primary 
vertex. It also excluded duplicate clusters resulting from the overlap of the pixel 
modules in <fr and from a small set of pixels at the centres of the pixel modules 
that share readout channels [12) . 

The high charged particle multiplicity in Pb+Pb collisions can generate 
misidentified tracks and/or two-point tracklets when only two or three mea- 
surements are made on each trajectory. The misidentified contributions have 
been evaluated using the MC studies described below, but to check the MC re- 
sults, an independent, data-driven estimate of misidentified two-point tracklets 
was obtained using a variant of the two-point tracklet algorithm. In the default 
two-point tracklet analysis, referred to as "Method 1" , at most one tracklet was 
reconstructed for a given cluster on the first pixel layer. If multiple clusters 
on the second pixel layer fell within the search region defined in Equation [T] 
the closest cluster to the projected position was chosen. This method limits, 
but does not eliminate, the generation of misidentified tracklets. A second im- 
plementation of the two-point tracklet algorithm, referred to as "Method 2", 
produced tracklets for all combinations of clusters on the two layers consistent 
with the search region. Using Method 2, the rate of false tracklets resulting from 
random combinations of clusters was estimated by performing the same analy- 
sis but with the clusters on the second layer having their z positions inverted 
around the primary vertex and their azimuthal angles inverted, (f> — > n — <f>. The 
tracklet yield from this "flipped" analysis was then subtracted from the proper 
tracklet yield event-by-event to obtain the estimated yield of true tracklets, 



where iVfp represents the yield of two-point tracklets using Method 2 and N% 
represents the yield obtained by flipping the clusters in the second pixel layer. 
For the 0-10% centrality interval, the flipped yield is about 50% of the unflippcd 
yield in the I77I < 0.5 region. 

The response of the detector to the charged particles produced in Pb+Pb 
collisions and the performance of the track and tracklet methods was evaluated 
by MC simulations of Pb+Pb collisions using the HIJING [5] event generator 
followed by GEANT4 [TO] simulations of the detector response [2D] . The result- 
ing events were then reconstructed and analyzed using the full offline analysis 
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chain that was applied to the experimental data. HIJING events were gener- 
ated without jet quenching and with an unbiased impact parameter distribution. 
Impact parameter and px-dependent elliptic flow was imposed on the HIJING 
events after generation and prior to simulation. The GEANT4 detector ge- 
ometry included a distribution of disabled pixel modules matching that in the 
experiment. The MC events were used to derive correction factors from recon- 
structed pixel tracks and two-point tracklets to the primary HIJING particles. 
Primary particles were defined to be either particles originating directly from 
the Pb+Pb collision or particles resulting from secondary decays of HIJING 
produced particles with lifetimes cr < 1 cm. 

From the MC simulated events, correction factors accounting for particle 
detection efficiency, misidcntified tracks or tracklets from unrelated clusters, 
and extra tracks or tracklets from secondary decays or from interactions in the 
detector were calculated. The correction factors were evaluated in 20 intervals 
of detector occupancy (O) parameterized using the number of reconstructed 
clusters in the first pixel layer in the region < 1. Different corrections were 
applied to the pixel track and both two-point tracklet measurements. For the 
pixel tracks, the efficiency, e pt , for reconstructing tracks associated with charged 
primary particles was obtained from 



where N pr represents the number of charged primary particles produced by 
HIJING within a given 77 interval, and N™ atch represents the portion of those 
primary particles matched to reconstructed pixel tracks. The contributions to 
the number of reconstructed pixel tracks (AT pt ) from "background" sources were 
separately evaluated to produce a "background" fraction 

b ^°^= N pt (0,r,) ' < 4 > 

where 7V pt acks represents the number of tracklets from secondary interactions 
and decays, from particles initially produced outside the kinematic acceptance 
of the measurement but scattering into it, and from combinations of clusters not 
associated with any primary or secondary particle in the GEANT4 simulation. 
This factor was combined with e pt (0,rj) to produce a correction factor 

C pt (0, V )= £ J ^ (l-b pt (0, V )). (5) 

For the 0-10% centrality interval, e pt is about 0.55 and 6 pt is about 0.02 in the 
mid-rapidity region, giving a C p t of about 1.8. 

For the two-point tracklet methods, a single multiplicative correction factor 
was obtained from the MC simulations, 

r (O n) - MglZ?) (ft) 

C2 ^ v) = N^fy (6) 
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Figure 1: Tracklet candidate Ar] (left) and A</> (right) distributions from data (histogram) 
and reweighted MC (shaded region) for Pb+Pb collisions at v / s^7 = 2.76 TeV. The top 
panels correspond to |r?| < 1 and the bottom panels correspond to 1 < \rj\ < 2. Data and MC 
distributions are normalized to the same area. 
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where N 2p (0,ri) represents reconstructed tracklcts. For the two-point tracklet 
Method 2, N2 P (0, rf) was obtained from the MC events via Eq. fusing the same 
flipping procedure as that applied in the data. For the 0-10% centrality interval, 
the correction factor is about 1.05 for Method 1 and 1.25 for Method 2 in the 
mid-rapidity region. 

The Pb+Pb charged particle px spectrum measured at ^/s NN = 2.76 TeV 
[2"T] differs from the spectrum generated by HIJING at low and high px, with 
the generator exceeding the data by 20% at px = 500 MeV, and underpre- 
dicting the charged particle yield by a factor of about two at px = 1-5 GeV. 
Because the MC corrections are applied to the data in matching O intervals, the 
mismatch in the spectrum does not influence the corrections for misidcntificd 
tracks or occupancy-induced inefficiencies. However, if left uncorrected the mis- 
match could distort the px-weighted single track or tracklet efficiencies in the 
calculated correction factors. To avoid this distortion a px-dependent weight 
was applied to the generated particles and to tracklets or tracks that match 
generated particles in Eqs. |3][6j The px-dependent weights were obtained us- 
ing an iterative procedure that, in each analyzed centrality interval, optimally 
matched the px spectrum of pixel tracks in Pb+Pb data with the solenoid mag- 
net turned-on to the reweighted spectrum produced from a separate sample 
of HIJING+GEANT4 simulations also performed with the solenoid turned-on. 
Distributions of Arj and A</> for candidate tracklets are shown in Fig. [T] for two 
different pseudorapidity intervals, \t)\ < 1 and 1 < \r}\ < 2. The corresponding 
distributions for the reweighted HIJING+GEANT4 events are also shown in 
the figure and compare well with the data. The maximum difference between 
data and MC is less than 5%. It should be noted that the cr^rj) and 00(77) 
mentioned above are evaluated using the unreweighted MC, but they are ap- 
plied consistently to data and reweighted MC when calculating all ^-dependent 
corrections. 

Uncorrected pixel track and two-point tracklet pseudorapidity distributions 
for 0-10% centrality collisions are shown in the top left panel of Fig. [2] The 
corrections described above are applied to obtain corrected, per-event primary 
charged particle pseudorapidity distributions, averaged over the events in each 
centrality bin (c), according to 



dN ch 
di] 



1 A M raw 



where AiV raw indicates either the number of reconstructed pixel tracklets or 
two-point tracklets and C(0,rf) indicates the 77-dependent correction factors 
corresponding to the occupancy bin for each event. The corrected dN^/dr] dis- 
tributions for the 0-10% centrality interval are shown in the middle left panel of 
Fig. [2] The bottom left panel of Fig. [2] shows the ratio of the pixel tracking and 
two-point tracklet Method 2 results to the two-point tracklet Method 1 results. 
In spite of the factor of ^2 differences between the raw yields for the three 
reconstruction methods, the corrected pseudorapidity distributions for central 
collisions agree within 5%. The measurements presented in the remainder of 
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Figure 2: Left: Top: uncorrected track/tracklet d/V raw /iir/ distribution from tracklet 
Method 1 (points), tracklet Method 2 (squares) and pixel tracking (blue triangles) for 0- 
10% centrality events. Middle: corrected tracklet and track dN^/d-q distributions. Bottom: 
ratio of dN c ^/dr) from the tracklet Method 2 (squares) and pixel tracking (triangles) to track- 
let Method 1. Right: dN^/dr) distributions from tracklet Method 1 for eight 10% centrality 
intervals. The statistical errors are shown as bars and the systematic errors are shown as 
shaded bands. 
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this paper were obtained from tracklet Method 1, which has the highest recon- 
struction efficiency, only a moderate contribution of misidentihed tracklets, and 
the smallest correction factors. The resulting corrected dN^/dr] distributions 
are shown for 8 centrality intervals in the right-hand panel of Fig. [2] 

5. Systematic uncertainties 

Various studies were performed to quantify the experimental uncertainties 
on the dNch/dr) measurement. To address inaccuracies in the MC description of 
bad channels, disabled sensors, or other small instrumental problems, a compar- 
ison was made of unit-normalized rj and <f> distributions of clusters in each of the 
first two pixel layers between data and MC. The agreement between the rj and 
4> distributions was found to be better than 0.05% and 0.4% in the first and sec- 
ond layers, respectively. Therefore, a combined systematic uncertainty of 0.4% 
is assigned to account for potential MC inaccuracies. To evaluate the impact of 
inaccuracies in the description of the detector material in the GEANT4 simu- 
lation, a separate set of HIJING+GEANT4 simulations was performed with an 
artificial 10% increase in detector material and a 15-20% increase in material in 
various non-instrumented regions. The results obtained using correction factors 
from this "extra material" sample agree with those obtained using the default 
corrections to better than 2%. Furthermore, the analysis was repeated using 
a different ATZ selection (see Eq. [IJ, ATZ < 1.5, which should have a different 
sensitivity to multiple scattering, secondaries, and occupancy effects. The cor- 
rections for the AT? < 1.5 selection differ from those of the default analysis in 
central (0-10%) collisions by 10% at rj = and 20% at rj = 2. However, the cor- 
rected pseudorapidity distributions agree to 1% in all centrality intervals. To 
address differences between the HIJING description of particle production in 
Pb+Pb collisions and reality, the analysis was performed without the px spec- 
trum re-weighting; the results agree with those obtained using the re-weighting 
within 0.5%. To address potential errors resulting from discrepancies in parti- 
cle composition between data and MC, the changes in correction factors that 
would result from enhanced charged kaon and proton production as observed 
at RHIC [55] have been evaluated. From the impact of the modified corrections 
on the final result, a 1% systematic uncertainty due to incomplete knowledge 
of the hadron composition is assigned. To further test the sensitivity of the 
results to the use of the HIJING generator, a set of MC simulations using the 
HYDJET event generator [23] was produced, and a separate set of correction 
factors was obtained from this MC sample. HYDJET has a more complete 
description of soft particle production than HIJING, including a description of 
elliptic flow, and the version used here was tuned to have much lower multi- 
plicities than found in HIJING. In central collisions, the results obtained using 
the HYDJET-based corrections agree with the HIJING-based results to bet- 
ter than 0.5% at mid-rapidity, but differ by as much as 7.5% at T] = ±2. A 
centrality-dependent and ^-dependent systematic error is assigned to account 
for this difference. To address the inaccuracies from the analysis procedure, 
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Source 


Uncertainty (0-10%) 


(70-80%) 


MC detector description 


0.4% 


0.4% 


Extra material 


2% 


2% 


ATI cut 


1% 


1% 


px re-weighting 


0.5% 


0.5% 


Hadron composition 


1% 


1% 


Enhanced K a , A 


1% 


1% 


HYDJET 


0.5-7.5% vs. r) 


0% 


Analysis Method 


3.5% 


1% 


Combined (r/ = 0) 


4% 


3% 


Combined (r; = 2) 


8.5% 


3% 



Table 1: Summary of the various sources of systematic uncertainties and their estimated 
impact on the dN c ^/drj measurement in central (0-10%) and peripheral (70-80%) Pb+Pb 
collisions. Only the uncertainty due to the choice of the event generator is jj-dependent. 



a systematic uncertainty is assigned based on the differences between the re- 
sults obtained from the three reconstruction methods described in this paper. 
That uncertainty is centrality-dependent and maximal for the 0-10% centrality 
interval for which a 3.5% uncertainty on the overall scale of the pseudorapid- 
ity distribution is assigned based on the comparison of the three results in the 
left, bottom panel of Fig. [2] The systematic uncertainties described above are 
summarized in Table [T] for the most central (0-10%) and the most peripheral 
(70-80%) intervals. The total systematic uncertainties are shown as shaded 
bands in the right panel of Fig. [2] 



6. Results 

The measured charged particle dN c h/dr] shown in Fig. [2j increases rapidly 
with collision centrality for all 77. It is conventional to characterize par- 
ticle production in nucleus-nucleus collisions by the mid-rapidity dN^/drj, 
dN c b/ 'd^—o, which here is defined to be dN c h/di] averaged over I77I < 0.5. 
The analysis presented in this paper yields dN^/ 'dr)\ v =o values in central 
collisions of 1479 ± 10(stat.) ± 63(syst.), 1598 ± ll(stat.) ± 68(syst.), and 
1738 ± 12(stat.) ± 75(syst.) for the 0-10%, 0-6%, and 0-2% centrality inter- 
vals, respectively. Table [2] provides results of the dN c ^/ 'rfr/l T?=0 measurements 
for all centrality bins. 

The top panel of Fig. [3] compares the ATLAS measurement to the pre- 
viously reported ALICE [S] and CMS results for \r)\ < 0.5 for the 0- 
5% centrality interval in terms of dN c h/dr]\ v= o per colliding nucleon pair, 
dA 7 c h/d?7|r ;= o/(( Ap ar t )/2), and to other A+A measurements at different y/s NN 
(see 0, which includes data from Refs.[H]-[2!]). The ALICE and CMS 0-5% 
centrality measurements agree with the result reported here for the 0-6% cen- 
trality interval, 8.5 ± 0.1(stat.) ± 0.4(syst.), within the quoted errors. The LHC 
results show that the multiplicity in central A+A collisions rises rapidly with 
^/s NN above the RHIC top energy of ,/s NN =200 GeV. The three curves shown in 
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Ccntrali ty 


< Np ar t > 


<i7V ch /d»?|i7=o 


dN ch /dri\ v=0 /{N pa , It }/2 


0-2% 


396 ± 2 


1738 ± 76 


8.8 ± 0.4 


2-4% 


378 ± 2 


1591 ± 67 


8.4 ± 0.4 


4-6% 


356 ± 3 


1467 ± 63 


8.2 ± 0.4 


6-8% 


335 ± 3 


1350 ± 57 


8.1 ± 0.4 


8-10% 


315 ± 3 


1250 ± 53 


8.0 ± 0.3 


10-12% 


296 ± 3 


1159 ± 48 


7.8 ± 0.3 


12-14% 


277 ± 4 


1074 ± 44 


7.8 ± 0.3 


14-16% 


260 ± 4 


996 ± 41 


7.7 ± 0.3 


16-18% 


243 ± 4 


918 ± 37 


7.6 ± 0.3 


18-2070 


228 ± 4 


849 ± 34 


7.5 ± 0.3 


20-25% 


203 ± 4 


739 ± 29 


7.3 ± 0.3 


25-30% 


170 ± 4 


603 ± 24 


7.1 ± 0.3 


30-35% 


142 ± 4 


486 ± 19 


6.9 ± 0.3 


35-40% 


117 ± 4 


387 ± 15 


6.6 ± 0.3 


40-45% 


95.0 ± 3.7 


303 ± 11 


6.4 ± 0.3 


45-50% 


76.1 ± 3.5 


233 ± 9 


6.1 ± 0.4 


50-55% 


59.9 ± 3.3 


176 ± 6 


5.9 ± 0.4 


55-60% 


46.1 ± 3.0 


129 ± 5 


5.7 ± 0.4 


60-65% 


34.7 ± 2.7 


93 ± 3 


5.3 ± 0.5 


65-70% 


25.4 ± 2.3 


65 ± 2 


5.1 ± 0.5 


70-75% 


18.0 ± 2.0 


43 ± 2 


4.8 ± 0.6 


75-80% 


12.3 ± 1.6 


28 ± 1 


4.6 ± 0.6 



Table 2: Tabulation of measurements of dN c ^/ dri^—Q evaluated over \r]\ < 0.5 and 
dN c h/dr)\ri=o/((Npaxt)/2) for the full set of centrality bins considered in the analysis and 
shown in Fig. [3] The uncertainties on dN c ^/drj\ v= o include statistical and systematic er- 
rors on the multiplicity measurement. The errors reported for dN c -^/ dr}^— o/((JVpart)/2) also 
include systematic uncertainties on the centrality selection and (jVp ar t) determination. 



Fig. [3 indicate possible variations of dN c i 1 /dr]\ 71=0 /((Np aT t} /2) with ^/s NN . The 
dotted curve describes a ^/s NN dependence expected from Landau hydrodynam- 
ics [7] . It is clearly inconsistent with the data. The dot-dashed curve represents 
a logarithmic extrapolation of RHIC and SPS data [30 that is also excluded by 
the measurement presented in this paper and by the ALICE and CMS measure- 
ments. The dashed curve shows an s 015 dependence suggested by ALICE [8 
that is consistent with the ATLAS measurement. Also shown in the top panel 
in Fig. [3] are results from p+p and p+p measurements at different y/s ([2] and 
references therein, as well as [3T]-[3S]). The excess of dN^/ dri\ v=0 / {{N vart ) /2) 
in A+A collisions over p+p collisions observed at RHIC persists and is propor- 
tionately larger at the higher ^/s NN values of the LHC. 

The bottom panel of Fig. [3] shows dN c h,/dr]\ ri —o/({Np ax t)/2) as a func- 
tion of (iVp art ) for 2% centrality intervals over 0-20%, and 5% centrality in- 
tervals over 20-80%. The values are also reported in Table [2j A mod- 
erate variation of dN c t 1 /dri\ ri —o/((N paI t)/2) with (A part ) is observed, from a 
value of 4.6 ± 0.1(stat.) ± 0.6(syst.) at (Apart) = 12.3 (centrality 75-80%) to 
8.8 ± 0.1(stat.) ± 0.4(syst.) at (A part ) = 396 (centrality 0-2%). The increase 
of dN c h/dT)\ v= o/((N piiIt )/2) with (A par t) is monotonic up to the most central 
interval (0-2%). This demonstrates that, even for the most central collisions, 
variations in centrality - as characterized by transverse energy depositions well 
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Figure 3: Top: ^/s NN dependence of the charged particle dN c ^/dr] per colliding nucleon 
pair dN c i L /dr]\n=o/({Npa,Tt)/i) from a variety of measurements in p+p and p+p (inelastic and 
non-single diffractive results from 2 and references therein, as well as 31 - 35 ) and central 
A+A collisions, including the ATLAS 0-6% centrality measurement reported here for < 0.5 
and the previous 0-5% centrality ALICE [8] and CMS [9] measurements (points shifted hor- 
izontally for clarity). The curves show different expectations for the y^ s nn dependence in 
A+A collisions: results of a Landau hydrodynamics calculation (dotted line) , an s 015 
extrapolation of RHIC and SPS data proposed by ALICE [8] (dashed line), a logarithmic ex- 
trapolation of RHIC and SPS data from [30] (solid line). Bottom: dN ch /dr]\ v= o/{(Np a , It )/2) 
vs (JVpart) for 2% centrality intervals over 0-20% and 5% centrality intervals over 20-80%. 
Error bars represent combined statistical and systematic uncertainties on the dN^/dr]]^— q 
measurements, whereas the shaded band indicates the total systematic uncertainty including 
(^Vpart) uncertainties. The RHIC measurements (see text) have been multiplied by 2.15 to 
allow comparison with the y/s NN = 2.76 TeV results. The inset shows the (AT par t) < 60 region 
in more detail. 



13 



outside the acceptance used for the multiplicity measurement - yield significant 
changes in the measured final state multiplicity. 

The bottom panel of Fig. [3] also shows ALICE and CMS measurements 
of dNch/drjlrj^o as a function of (iV part ) that agree with the results presented 
here for all centrality intervals. Also shown are results from Au+Au collisions at 
^/s NN = 200 GeV obtained from an average of measurements from the four RHIC 
collaborations j3S]-[in]. Similar to the approach used in Ref. [5], the 200 GeV 
Au+Au results have been scaled by a factor of 2.15 to allow comparison with the 
^/s NN =2.76 TeV data. This factor was obtained by matching the most central 
200 GeV Au+Au dN ch /dr] measurement at r] = to the dN ch /dr] measurement 
from this paper at 77 = in the 2-4% centrality interval, the interval that has 
the closest value of (V part ) to the most central 200 GeV measurement. After re- 
scaling, the trend of the 200 GeV data is in good agreement with the 2.76 TeV 
measurements for all reported centrality intervals. Similar observations have 
been made previously in comparisons of top energy RHIC data to much lower 
energies [2]. Therefore, this scaling behavior appears to be a robust feature of 
particle production in heavy ion collisions. 

To evaluate the shapes of the measured charged particle dN^/drj distribu- 
tions Fig. Pfl (top) shows the dN^/drj distribution divided by dN^/ drj\ n= Q for 
the 70-80% centrality interval. For this centrality interval, the dN^/drj in- 
creases by 7% ± 1% from r\ = to |r;| > 1. The bottom panel shows ratios 
of d V c h/d?7/(( V part )/2) for several other 10% centrality intervals to the same 
quantity in the 70-80% interval. No significant variation of the shape of dN^/dr] 
with centrality is observed within the systematic uncertainties. 



7. Conclusions 

This paper presents results on the measurement of charged particle pseudo- 
rapidity distributions over |7y| < 2 as a function of collision centrality in a sample 
of ^/s NN =2.76 TeV lead-lead collisions recorded with the ATLAS detector at 
the LHC. Three different analysis methods are used, based on the pixel detec- 
tor and using events with the solenoid magnet turned off in order to measure 
particles with transverse momenta as low as 30 MeV. The charged particle mid- 
rapidity dN^/dr), normalized by (V part )/2, is found to increase significantly 
with beam energy by about a factor of two relative to earlier RHIC data, and 
is substantially larger than p+p data at the same energy. The relative central- 
ity dependence of dN c h/drj\, 1 =o/((N par t) /2) agrees well with that observed at 
RHIC. These results agree well with previous mid-rapidity measurements from 
ALICE and CMS. Furthermore, the peripheral (70-80%) dN ch /drj distribution 
shows a significant rise with increasing \r)\ away from r\ = 0. No variation 
of the shape of the dN^/drj distribution with centrality outside the reported 
systematic uncertainties is observed. 
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Figure 4: Top: dN c ^/dr] distributions from tracklet Method 1, scaled by dN c h/ 'dr]\ v —o, as a 
function of the pseudorapidity for the 70-80% centrality interval. The statistical errors are 
shown as error bars. Bottom: Ratio of dN c ^/ dr) / ({-/V par t} /2) measured in different central- 
ity intervals: 0-10% (squares), 20-30% (triangles), 40-50% (inverted triangles) and 60-70% 
(crosses) to that measured in peripheral collisions (70-80%). Statistical uncertainties are 
shown as bars while jj-dependent systematic uncertainties are shown as shaded bands. 
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